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We extend a recently introduced method for computing Casimir forces between arbitrarily-shaped 
metallic objects [M. T. H. Reid et ai, Phys. Rev. Lett. 103 040401 (2009)] to allow treatment of 
objects with arbitrary material properties, including imperfect conductors, dielectrics, and magnetic 
materials. Our original method considered electric currents on the surfaces of the interacting ob- 
jects; the extended method considers both electric and magnetic surface current distributions, and 
obtains the Casimir energy of a configuration of objects in terms of the interactions of these effective 
surface currents. Using this new technique, we present the first predictions of Casimir interactions in 
several experimentally relevant geometries that would be difficult to treat with any existing method. 
In particular, we investigate Casimir interactions between dielectric nanodisks embedded in a di- 
electric fluid; we identify the threshold surface-surface separation at which finite-size effects become 
relevant, and we map the rotational energy landscape of bound nanoparticle diclusters. 

PACS numbers: 03.70.+k, 12.20.-m, 42.50.Lc, 03.65.Db 



Rapid progress in experimental Casimir physics pQ has 
created a demand for theoretical tools that can predict 
Casimir forces and torques on bodies with realistic ma- 
terial properties, configured in experimentally relevant 
geometries. Techniques for computing Casimir interac- 
tions have traditionally fallen into two categories [2J. In 
the scattering approach [3], the Casimir energy is ob- 
tained from spectral properties of the classical electro- 
magnetic (EM) scattering matrix of the geometry; this 
method leads to rapidly convergent and even analyti- 
cally tractable 0] formulae for certain highly symmet- 
ric geometries in which a special choice of coordinates 
leads to a simple form for the scattering matrix, but is 
less readily applicable to objects of general asymmetric 
shapes. In the alternative numerical stress-tensor (ST) 
approach [5] , the Casimir force and torque on a compact 
body are obtained by numerically integrating the (ther- 
mally and quantum-mechanically averaged) Maxwell ST 
over a bounding surface surrounding that body, with the 
value of the ST at each spatial point obtained by solving 
a classical EM scattering problem. The ST method has 
the virtue of flexibility, as almost any of the myriad ex- 
isting methods for computational EM packages may be 
applied off-the-shelf to solve scattering problems for ar- 
bitrary geometries and materials, but the ST integration 
adds a layer of conceptual and computational complexity. 

In recent work [6], we have introduced a new method 
for Casimir computations that combines the best fea- 
tures of the scattering and stress-tensor approaches. Our 
surface-current technique discretizes the surfaces of in- 
teracting objects into small surface patches (Fig. 1) — 
thus retaining the full flexibility of the numerical ST 
approach in handling arbitrarily complex, asymmetric 



geometries — but bypasses the unwieldy ST integration 
by obtaining the Casimir energy directly from spec- 
tral properties of a matrix describing the interactions 
of surface currents on the discretized object surfaces — 
thus mirroring the conceptual directness of the scatter- 
ing approach, dramatically improving computational effi- 
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FIG. 1. Casimir force between silicon and teflon spheres and 
cubes immersed in ethanol, normalized by the PFA predic- 
tions (see text). For the sphere-sphere case, the solid circles 
represent data computed using the methods described in this 
paper, while the thick solid line reproduces data from Ref. [?]• 
For the other three cases, the solid and hollow circles represent 
data computed using the methods described in this paper, and 
the dotted lines are a guide to the eye. 
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ciency for moderate-size problems where dense linear al- 
gebra is applicable, and establishing a connection to the 
equivalent-surface-current method of classical EM [8] . 

The technique discussed in [5] was limited to the case of 
perfect-metal bodies in vacuum. In this paper we extend 
the method to allow treatment of materials with arbitrary 
frequency-dependent (but spatially constant) permittiv- 
ity e and permeability fi. A thorough description of our 
extended formulation will be provided elsewhere [5] ; here 
we present a brief overview of our method and demon- 
strate its use in predicting Casimir interactions in sys- 
tems that would be difficult to treat by any other method. 

Casimir Interactions from Surface Currents. The 
Casimir energy of a collection of objects is given quite 
generally by an expression of the form [TU] 



h 

'27r. /M 



d£ log 



(1) 



where the partition function for the electromagnetic field 
at imaginary frequency £ is 

= J [DA] e ~i / A (5' x )-5" 1 (f.x,x')-A(f,x')<ixdx'. 

here Q is the (gauge-fixed) photon propagator, and the 
notation [• • -] c indicates that this is a constrained path 
integral, in which the path integration ranges only over 
field configurations A that satisfy the boundary condi- 
tions in the presence of the interacting objects. {Z^ in 
is Z with all objects removed to infinite separations.) 
Although the exponent in ([2| is only quadratic in the 
integration variables, the implicit constraint in the in- 
tegration measure prevents direct evaluation of the path 
integral. We obtain a more tractable expression by repre- 
senting the constraints on the A field explicitly via func- 
tional 5 functions [TTJ [T2] . Consider a single point x on 
the surface of an object in our geometry. The boundary 
conditions at x are that the components of E and H tan- 
gential to the object surface are continuous as we pass 
from the interior to the exterior of the object at x: 

Ejp (x) = Ejj ut (x) , H}| n (x) = Hff ut (x) . 

We introduce two 2-dimensional S functions that impose 
these constraints at x: 
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where we may think of the Lagrange multipliers K and N 
as two-component vectors living in the tangent space to 
the object surface at x. Aggregating the corresponding 
5 functions for all points x on all object surfaces, we 
obtain an explicit functional S function constraining the 
electromagnetic field to satisfy the boundary conditions: 



A] =/^ K / PNe^ lK£EA+lN£HA }' 
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where the integration in the exponent is over the sur- 
faces of all objects in our geometry, the path integrations 
range over all possible tangential vector fields K(x), N(x) 
on the object surfaces, and £ E ,£ H are differential oper- 
ators that act on A to yield E and H. Because K and 
N respectively enforce the boundary conditions on the 
electric and magnetic fields, it is tempting to interpret 
these tangential vector fields as electric and magnetic 
surface current distributions on the surfaces of the in- 
teracting objects. Effective surface current distributions 
of this sort are ubiquitous in equivalence-principle tech- 
niques for EM field computation [Sj and form the basis of 
the boundary-element method of computational EM [13] . 

Inserting the explicit constraint ([3| into the path inte- 
gral ^ allows us to relax the implicit constraint [• • • ] c , 
whereupon the path integral over the A field may be 
performed analytically: 
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where ~ denotes equality up to constant factors that can- 
cel in the ratio in ([I]), and where the effective action 
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describes the interactions of electric and magnetic surface 
currents on the bodies of the interacting objects. The 
frequency-dependent material parameters e(£), /x(£) of all 
interacting objects, and of the medium in which they are 
embedded, enter into through £ E ,£ H , and £/. 

Discretization. To evaluate Q for a given collection 
of objects, we now proceed as in [5] by discretizing the 
surfaces of the objects into small planar triangles [14] . 
We introduce a set [15] of localized vector-valued ba- 
sis functions {f a } and approximate the surface currents 
K(x),N(x) as expansions in the finite set {f a }: 

K(x)«^if Q f a (x), N(x)«^A^ Q f Q (x). 

a a 

The functional integral in Q becomes a finite- 
dimensional Gaussian integral, 

J PKPNe^M- j J[dK a dN a e-O m ^\%) 

~ [detM(0]r 1/a (5) 
and, inserting into 0, the Casimir energy becomes 
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while differentiating with respect to the position of an 
object gives the Casimir force on that object: 
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The matrix M describes interactions among the effective 
electric and magnetic surface currents. More specifically, 
each pair of localized basis functions (f a , fg) corresponds 
to a 2x2 block of M of the form 



M a/3 = 



where e.g. G ME is a dyadic Green's function representing 
the magnetic field due to a point electric dipole source. 

Comparison to the BEM. The matrix M has precisely 
the same form as the matrix that arises in the boundary- 
element method (BEM) for electromagnetic scattering 
problems |13| . Indeed, once we have implemented the 
computational machinery needed to evaluate ([7]), we 
could alternatively use this same machinery to compute 
Casimir forces in a different way — namely, by numeri- 
cally integrating the ST over a bounding surface, with 
the value of the ST at each spatial quadrature point ob- 
tained by solving a BEM scattering problem [16 . This 
procedure would entail solving linear systems M • X = B 
for multiple vectors B, which in practice would be han- 
dled by first decomposing M into triangular factors. Such 
a factorization step is also required for numerical evalu- 
ation of the determinant or trace expressions in (6|7|. 
However, in our formulation, this factorization step is 
all that is required, with the Casimir energy and force 
then obtained immediately from the simple relations ^ 
[7]); in contrast, the ST method requires extensive pre- 
and post-processing, in addition to the factorization of 
the BEM matrix, to evaluate the ST integral. Thus, in 
addition to the conceptual simplicity of our equations 
([6] [7j) as compared to the ST formulation, our method 
is significantly more computationally efficient than the 
ST approach, at least as long as we consider problems 
of moderate size (N = dim M < 10 4 ), for which direct 
linear-algebra methods are feasible. 

As an explicit demonstration of this computational ad- 
vantage, the three points represented by hollow circles on 
the cube-cube data curve in Figure 1 were computed by 
the ST-integration method using our BEM code. The 
resulting force predictions agree with the predictions of 
Q, but consume some 12-24 times more CPU time. 

Moreover, although here we derived our equation Q in 
the context of the scattering approach to Casimir physics, 
the same equation can alternatively be shown to follow 
directly from an application of the BEM to stress-tensor 
Casimir computations, with the spatial integral evalu- 
ated analytically. Full details of these two complemen- 
tary derivations will be presented in [5]. 

Applications. We now use our extended method to 
predict Casimir interactions in a number of new systems 
that would be difficult to treat with any other method. 

Attractive and repulsive Casimir forces between sili- 
con and teflon spheres and cubes immersed in ethanola. 
Ref. [7J , using the methods of Ref . [I] , predicted Casimir 



forces between spheres and slabs composed of various ma- 
terials and suspended in ethanol. Figure 1 validates our 
method by using it to reproduce the results of Ref. [7J for 
the Casimir force between teflon and (undoped) silicon 
spheres, then extends these results by plotting Casimir 
forces between spheres and cubes of the same material 
properties (with the cubes sized to have the same vol- 
ume as the corresponding spheres). For each geometry, 
we plot the ratio of the computed Casimir force to the 
force predicted by the perfect-metal proximity force ap- 
proximation (PFA) for that geometry; since the PFA al- 
ways predicts an attractive force, a positive (negative) 
ratio implies an attractive (repulsive) Casimir force. It is 
interesting that all four geometries exhibit a repulsive- 
attractive force crossover at approximately the same 
surface-surface separation distance. 

Finite-size effects in nanoparticle diclusters. Ref. [7] 
predicted Casimir forces between dielectric slabs, of infi- 
nite cross-sectional-area, embedded in ethanol. It is in- 
teresting to ask how these predictions are modified for 
finite slabs (discs). Figure 2 plots the PFA- normalized 
Casimir force between pairs of puck-shaped nanoparti- 
cles, one composed of polystyrene and the other of doped 
silicon, immersed in ethanol. For comparison, the dotted 
curves indicate the (PFA-normalized) force between pairs 
of infinite slabs of the same materials. For small values 
of the surface-surface separation d, the normalized puck- 
puck force tends to the slab-slab limit, but finite-size ef- 
fects cause the curves to deviate beyond separations of 
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FIG. 2. Casimir force between puck-shaped polystyrene and 
doped-silicon (dopant density p) nanoparticles immersed in 
ethanol. The dotted curves indicate the force on slabs of 
the same materials and thickness but infinite cross-sectional 
area [7]. (Material properties are taken from Ref. [7].) 
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FIG. 3. (Color) Rotational stability of hockey-puck shaped nanoparticles suspended in ethanol. One nanoparticle is composed 
of polystyrene, and the other of silicon with a dopant concentration of 1.4- 10 19 cm -3 (corresponding to the middle curve in the 
previous plot). We fix the centers of the particles at center-center separations of S = {500, 700} nm allow them to rotate about 
their centers. The minimum-energy configurations are depicted in the insets at right. (All angles are measured in degrees.) 



d ~ R/2 (with R the puck radius). 

Rotational stability of dielectric nanoparticles. For 
the nanoparticle pair corresponding to the intermediate 
curve in Figure 2 (dopant density of 1.4 • 10 19 cm ~ 3 for 
the silicon nanoparticle), we now fix the centroids of the 
nanoparticles and plot the Casimir energy (Figure 3) as 
a function of rotation angles (see the inset of Figure 3) 
for two centroid-centroid separation distances S of 700 
and 500 nm. The slight shift in the value of S changes 
the stable equilibrium (right). 
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